  function [freqv,perv,edgew]=freq_linspace(lb,ub,npoints)
% ========================================================== 
wlb=2*pi/ub;wub=2*pi/lb; 
freqv=linspace(wlb,wub,npoints+2); 
freqv=freqv(2:end-1);
perv=(2*pi)./(freqv); 
edgew=[(freqv(1)-wlb) (wub-freqv(end))];
if any( edgew > ((freqv(2)-freqv(1))+1e-5)   );
    warning('Egde bound widths might be too large');
end 